Dimerization-assisted energy transport in light-harvesting complexes 
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We study the role of the dimer structure of Hght-harvesting complex II (LH2) in excitation transfer from the LH2 
(without a reaction center (RC)) to the LHl (surrounding the RC), or from the LH2 to another LH2. The excited and 
un-excited states of a bacteriochlorophyll (BChl) are modeled by a quasi- spin. In the framework of quantum open 
system theory, we represent the excitation transfer as the total leakage of the LH2 system and then calculate the transfer 
efficiency and average transfer time. For different initial states with various quantum superposition properties, we study 
how the dimerization of the B850 BChl ring can enhance the transfer efficiency and shorten the average transfer time. 



I. INTRODUCTION 

To face the present and forthcoming global energy crisis, 
human should search for clean and effective energy source. 
Recently the investigations on the basic energy science for 
this purpose has received great attention and experienced im- 
pressive progress based on the fundamental physics In 
photosynthetic process, the structural elegance and chemi- 
cal high efficiency of the natural system based on pigment 
molecules in transferring the energy of sunlight have stim- 
ulated a purpose driven investigation^"^ ^, finding artificial 
analogs of porphyrin-based chromophores. These artificial 
systems replicate the natural process of photosynthesis^ so 
that the much higher efficiencies could be gained than that 
obtained in the conventional solid systems^. It is because 
one of the most attractive features of photosynthesis is that 
the light energy can be captured and transported to the reac- 
tion center (RC) within about lOOps and with more than 95% 
efficiency"^' 

Actually, in most of the plants and bacterium, the pri- 
mary processes of photosynthesis are almost in common^^i^ii^: 
Light is harvested by antenna proteins containing many chro- 
mophores; then the electronic excitations are transferred to the 
RC sequentially, where photochemical reactions take place 
to convert the excitation energy into chemical energy. Most 
recent experiments have been able to exactly determine the 
time scales of various transfer processes by the ultra-fast laser 
technology These great progresses obviously offer us a 
chance to quantitatively make clear the underlying physical 
mechanism of the photosynthesis, so that people can con- 
struct the artificial photosynthesis devices in the future to 
reach the photon-energy and photon-electricity conversions 
with higher efficiency. For example, quantum interference ef- 
fects in energy transfer dynamics has been studied for the 
Fenna-Matthews-Olson (FMO) protein complex, and it was 
found^ that, for such molecular arrays, the spatial correla- 
tions in the phonon bath and its induced decoherence could 
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affect on the efficiency of the primary photosynthetic event. 
The present paper will similarly study the influences of spa- 
tial structure on the primary processes of photosynthesis for 
the light-harvesting complexes II (LH2). 

In the past, by making use of the x-ray crystallographic 
techniques, the structure of light-harvesting system has been 
elucidated^' In the purple photosynthetic bacteria, there ex- 
ist roughly two types of light-harvesting complexes, referred 
to as light-harvesting complex I (LHl) and light-harvesting 
complex II (LH2). In LHl, the RC is surrounded by a B875 
bacteriochlorophyll (BChl) ring with maximum absorption 
peak at 875 nm. The LH2 complex, however, does not con- 
tain the RC, but can transfer energy excitation to the RC in- 
directly through LHl. In the purple bacteria, LH2 is a ring- 
shaped aggregate built up by 8 (or 9) minimal units, where 
each unit consists of an ofyS-heterodimer, three BChls, and one 
carotenoid. The ayS-heterodimers, i.e., ^-apoproteins and JS- 
apoproteins constitute the skeleton of LH2, while the BChls 
are embedded in the scaffold to form a double-layered ring 
structure. The top ring including 16 (or 18) BChl molecules 
is named as B850 since it has the lowest-energy absorption 
maximum at 850 nm. The bottom ring with 8 BChls is called 
B800 because it mainly absorbs light at 800 nm. In every min- 
imal unit, the carotenoid connects B800 BChl with one of the 
two B850 BChls. Excitation is transferred from one pigment 
to the neighbor one through the Foster mechanism^, while the 
electron is spatially transferred via the Marcus mechanism^^. 
Generally, it is independent of the global geometry configura- 
tion of the system. 

In the present paper, we will study the energy transfer pro- 
cedure in LH2 by considering the structure dimerization of 
the B850 ring. It has been conjectured that the dimerized 
inter-pigment couplings can cause the energy gap to protect 
the collective excitations^^. Indeed, like the the Su-Schrieffer- 
Heeger model for the flexible poly acetylene chain^\ the 
dimerization of the spatial configuration with the Peierls dis- 
torted ground state will minimize the total energy for the 
phonon plus electron. As it is well known, this model exhibits 
a rich variety of nonlinear phenomena and topological excita- 
tions including the topological protection of the quantum state 
transfer^^. Similarly, we will show that, when the B850 ring 
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II. MODEL SETUP 

The simplified model of LH2 is shown in Fig. [T] All the 
bacteriochlorophylls (big and small green squares) are mod- 
eled by the two-level systems with excited state le^f), ground 

state 1^^^^), and energy level spacing Q^- The raising and low- 
ering quasi-spin operators of the jth two-level system on the 
[c] ring is expressed as 



■[c] 



[c] 



(1) 



where [c] = [a] ([b]) denotes the B800 (B850) BChl ring. Ap- 
proximately, all the couplings are supposed to be of XY type^. 
This simplification enjoys the main feature of excitation trans- 
fer. The Hamiltonians 



FIG. 1. (color online) The model setup of the light-harvesting com- 
plex II constructed by 8 unit cells. The couplings between the neigh- 
boring quasi-spins in the B850 ring is dimerized as J2(l + S) and 
^2(1 - gi denote the nearest couplings between B850 BChls and 
B800 BChls, while g2 denote the next nearest couplings between 
B850 BChls and B800 BChls. (a) Illustration of the whole system 
with g2 = 0. (b) Detailed drawing of three unit cells and their non- 
local couplings, (c) Legends. 



in LH2 is dimerized the excitation transfer efficiency may be 
enhanced to some extent. 

Based on the open quantum system theory, we simply 
model the excited and un-excited states of a BChl pigment as 
a quasispin. The excitation transfer is represented by the total 
leakage from a LH2. Using the master equation, we calculate 
the efficiency of excitation transfer and the average transfer 
time in low temperature for various initial states with diff'erent 
superposition properties. The results explicitly indicate that 
the dimerization of couplings indeed enhances the quantum 
transport efficiency and shortens the average transfer time. 

This paper is organized as follows. In Sec. II, a double-ring 
XY model with N unit cells is presented to simulate the LH2 
system. In Sec. Ill, the energy transfer process is described by 
the quantum master equation. The transfer efficiency t] (t) and 
the average transfer time r are introduced to characterize the 
dynamics of the system. In Sec. IV, we represent the master 
equation in the momentum space and show that only the (k, k)- 
blocks of the density matrix are relevant to energy transfer. In 
Sec. V, it is found that the transfer efficiency r] (t) and the aver- 
age transfer time r of an arbitrary initial state can be obtained 
through the channel decomposition. Some numerical analysis 
of T]^^'^^ (to) and r^^'^^ for all the ^-channels are presented in 
Sec. VI. They show that a suitable dimerization of the B850 
BChl ring can enhance the transfer efficiency and shorten the 
average transfer time. Conclusions are summarized at the end 
of the paper. In Appendix A, we provide an alternative way to 
deal with the energy leakage problem. In Appendix B, detail 
derivations of transforming the master equation from the real 
space to the ^- space are given. The approximate solution of 
j-[A,k] ^ = and k = ±71 channel is shown in Appendix C. 
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with = 8, describe the excitations of the B800 and B850 
BChl rings, respectively. In the B850 BChl ring, the parameter 
^ 9^ characterizes the dimerization due to the spatial defor- 
mation of the flexible B850 BChl ring in LH2. The coupling 
constants of Hiy are dimerized as /2 (1 + ^) and J 2 {I - since 
the intra-unit and inter-unit Mg-Mg distance between neigh- 
boring B850 BChls may be difl'erent. The non-local XY type 
interaction 



7=1 



7=1 



(4) 



is used to describe the interaction between the B800 and B850 
BChl rings. 

In the single excitation case, the quasi-spin can be repre- 
sented with a spinless fermion with the mapping 



(5) 



from the spin space Vs = C®^^ to the sub space Vp of the 
Fermion Fock space spanned by 

7) = 10) I O = A, 5, C; J = 1 , 2, . . . , A^} . (6) 

Hereafter, let us represent the site index as (0, 7), where j 
refers to a unit cell shown in Fig. [TJ and O = A, B,C to a. 
position type inside the unit cell. In the subscripts, the site 



index (O, j) is written as O j for simplicity. The vacuum state 
of the Fermion system |0) corresponds to the state that all the 
quasi- spins are in their ground states, 

2N 

Then the total Hamiltonian Hs = Ha -\- Ht, -\- Hat of LH2 is 
mapped into 

Hs=Yj [^aA]Aj + a, {b]Bj + C]Cj)] 

M 

N 

^ J2KI ^ 5)B]Cj ^ (1 - S)C]Bj^i] 

^g2A] [Bj^i + Bj.i + Cj^i + Cj-i) + H.c.} . (8) 

In the present work, no multi-fermion interactions are consid- 
ered for simplicity. 

On the other hand, we use the Holstein-Primakoff 
transformation^^ to map the quasi-spin into bosons. The ex- 
citations of the BChls can be described by quasi- spins with 
the total angular momentum S . Then D = A,B,C can be re- 
garded as the annihilation operators of bosons for the Fock 
space spanned by 

{(Z)jp|0) \D = A,B,CJ=h " ,N; 

nDj = OA,-" ,2S}. (9) 

For S > 1/2, one local bacteriochlorophyll has more than one 
excited states. In this case, higher order coherence could be 
included for further generalization. 

In the following, we focus on the single excitation case. 
Then the temperature should be suitable to ensure there is no 
higher order excited state. 

III. TRANSFER EFFICIENCY AND AVERAGE TRANSFER 
TIME VIA THE MASTER EQUATION 

Next we consider the energy transfer from an initial state 
p(0) = Y,PAj,Ai(0)\AJ}{AJ\, (10) 

which is a coherent superposition or a mixture of those local 
states |A, j) on the B800 ring. As time goes by, the initial state 
will evolves a state distributing around both the B800 and the 
B850 rings. Since there exists a difference of chemical poten- 
tial. AQ = Q.a - energy is transferred between the two 
rings during the time evolution. For an isolated LH2 system, 
such energy transfer is coherent, namely, the system oscillates 
between the B800 and the B850 rings . However, when a LH2 
is coupled to a heat reservoir with infinite degrees of freedom, 
irreversible energy transfer occurs. As illustrated in Fig. [51 in 
the real photosynthetic system, energy is transferred from one 



3 




FIG. 2. (color online) The first excited LH2 is treated as an open sys- 
tem while the other LHs are regarded as heat reservoirs. The energy 
transfer process is equivalent to the excitation leakage from the B850 
BChl ring of the LH2 system to the environment. 

LH2 to another LH2 or LHl through the B850 ring^^. There- 
fore, we regard the first excited LH2 as an open system, and 
the sum of others as the a heat reservoir. The energy transfer 
now can be manipulated as the energy leakage from the B850 
ring to the environment. 

In order to describe such a procedure that the excitations are 
finally transferred from the B850 ring to the heat reservoir, the 
Markovian master equation 

^=-i[Hs,p\+£(p) (11) 

in the Lindblad form is employed for determining the time- 
evolution of the density matrix. Here two kinds of loss pro- 
cesses, dissipation and dephasing, are considered as Lindblad 
terms 

X (p) = ^ [Xdiss,; (P) + Xdeph,; (p)] . (12) 

We suppose that each quasi-spin on the B850 ring is cou- 
pled to an independent heat reservoir-, which reflects the local 
modes of phonons and other local fluctuations. Then the dis- 
sipation from the jth unit cell is described as 

Xdissj (P) = rj (Ojpo] - \ {oJOy,p}), (13) 

where {•, •} denotes the anti-commutator. Here, the sink rate Y j 
at the 7th point may be site dependent. For the dynamics con- 
strained on the subsystem described by O operators, the last 
term of Xdissj (p) gives contribution -YjPoj,oj to dpoj,oj/dt, 
thus dissipation results in the reduction of the total popula- 
tion. Therefore, the dissipation term Eq. ([T3l) represents the 
incoherent transfer of energy into the environment. 
On the other hand, the dephasing term reads 

Xdeph.i (p) = r'jY, (o]ojpo]Oj - \ |oJo,-,p)). (14) 

0=B,C 
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Compared with the dissipation term, the dephasing one 
-^dephj (p) does not contribute to any time local change of the 
probability distribution, i.e., the derivative of the diagonal el- 
ements of the density matrix is irrelevant to this term. Thus 
the total population Yj0^a,b,c TjJ Poj,oj would be conserved if 
only the dephasing term were present. However, the dephas- 
ing process is also incoherent since it make the nondiagonal 
elements of the density matrix tend to zero. 

The above two contributions force the LH2 system to ulti- 
mately reach a steady state Psteady = |0) (0| = Pv,v, namely, in 
the long-time limit, all excitations are sinked away. The same 
steady state is obtained from Eq. ([TTt in the super-operator 
form 



M[p], 



(15) 



where [p] denotes the column vector defined by all matrix el- 
ements in some order, and the super-operator M is determined 
by 



M[p] = [-i[Hs,p\^£(p)l 



(16) 



In this sense the steady state is just the non-trivial eigenstate 
of M with vanishing eigen-energy. Usually, from det M = 0, 
the steady state can be found. 

However, we are interested in the system dynamics on a 
short timescale, i.e., how soon can the excitations be trans- 
ferred from one LH2 to the other light-harvesting complexes? 
To this end, the transfer efficiency t] (t) is defined as the popu- 
lation pv,v (t) of the vacuum state |0) at time t, 

ri(t)=Pv,v(t)- (17) 
The corresponding master equation ([TT1) 

'-^ = j]rjiO\{BJpB]^CjpC])\0} 

N 

= T,^jTjPoJ,oj (18) 

j=l O^B,C 

means that only the first term of Xdiss j (p) contributes to the 
time derivatives of pv,v (0- The transfer efficiency is given by 
the integral of the above formula^"-. 



^^^)= Z^^- Z Poj^ojiOdf. 

^0 j^l O^B,C 

The average transfer time r is further defined as^ 
1 ^ 

'^^ j^l 0=B,C 

= i ''Tj^jTj Poj.ojiOdf, 



(19) 



where usually 



7] = lim T](t) = I. 



(20) 
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FIG. 3. (color online) Configuration of the density matrix of the N = 
8 system in the subspace expanded by {|0> , |0, with O - A,B,C 
and k - (In/S) x 1, 2, • • • , 8. An initial state localized in the (^1, /^i)- 
block can be evolved to other (k, ^2 +^-^1) -blocks (black hollow dot- 
dash squares). Only the diagonal (^,^)-blocks (green solid squares) 
are related to the average transfer time. 



Therefore, an efficient energy transfer requires not only a per- 
fect transmission efficiency t] but also a short average time r. 

In Appendix A, we present an equivalent non-Hermitian 
Hamiltonian method, which can also be utilized to study the 
dynamics of the open system. 



IV. y^-SPACE REPRESENTATION OF THE MASTER 
EQUATION 

In this section we present the space representation of the 
above master equation, so that we can reduce the dynamics of 
time evolution in some invariant subspace. If all the dissipa- 
tion and dephasing rates are homogeneous on the B850 BChl 
ring, i.e., Tj = F and F^ = F', the whole system has trans- 
lational symmetry. For each unit cell containing three BChls 
shown in Fig. [H we introduce the Fourier transformation. 



Ol = J-ye-'^0] 



(22) 



for O = A,B, C. Then in the ^-space the Hamiltonian ([5]) is 
represented as //^ -Yuk^k with 

Hk = 2/1 cos kAlAk + {(gi -h 2g2 cos k) (a^^^^ A^q) 
+ J2 [a + ^) + (1 - ^) e-'^] BlCk + H.c.}. (23) 



Here k are chosen as discrete values 



2nl 



k= — , for / = 1,2,- •• ,A^. 
N 



(24) 
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In the subspace of the single excitation plus the vacuum 
with the basis 



N 

the general density matrix is decomposed into 

P = Pv,v + ^ ^kiM + ^ (Pv,k + 'Pk,v) . 

kiM k 

where 

Pv,v = Pv,v |0) (0| 

is the vacuum block while 

Pki,k2 = POkuO'k2 h) {0\ k2\ 

0,0'^A,B,C 



(25) 
(26) 

(27) 
(28) 



is called the (/:i, /:2)-block. For fixed and ^2, Poki,0'k2 form 
a matrix 



' pAkiMi PAkiMi PAkuCk2 ' 
pBkiMi PBkiMi PBki,Ck2 
PCkiMi PCkiMi PCki,Ck2 > 



(29) 



The ^- space representation of the density matrix is illustrated 
in Fig. [3] for the A/^ = 8 system. 

In the ^-space, the master equation (fTTI) is reduced to 



dt 



- -i(^kiPki,k2 - Pki,k2^k2) 



0=B,C 



N 



-i (F + r) (Oj 0,,p,,,, +p,,,,0j^0,,)| (30) 

for all the 

^ = -//^)^pii,v - ^ (F + r) y olOkPk,v, 

^ 0=B,C 

^ = iPv,kHk - ^ (F + r) ^ Pv,^o^o, 



O^B,C 



for all the ^, and 



(31) 



(32) 



O^B,C 



The details of the calculation are shown in Appendix B. 

We notice that the equations about ^k,v and pv,yt are com- 
pletely decoupled fromp^^,y^2 andpv,v. It follows from Eq. ([3Ql) 
that when no dephasing exists, i.e., F' = 0, the (^i,fe)-block 
^ki,k2 is decoupled with other p^'^^y^^ for (^^ ^2) ^ (^i^ fe)- Thus 
^ki,k2 only evolves in the (^1, fe)-block. However, when the 
dephasing is present (F' ^ 0), the term 

Z ^Z^I^^^^^-^^^l^^^^^-^^ 

O^B,C k 



actually induces the coupling between the (^1, /:2) -block and 
the (k, k2 -\- k - ki)-hlock. The initial py^i,y^2 evolves 
to 'pk,k2+k-ki as time goes by. A typical example of 
(/:, k2+k- ^i)-blocks are shown by the black hollow dot-dash 
squares in Fig. O The momentum diff'erence ki - is con- 
served during the evolution since 



k2- ki = {k2+k - ki) - k. 



(34) 



In addition, Eq. (1321) means that only the (k, ^)-blocks of the 
density matrix result in energy transfer, which are marked by 
the 8 green solid squares in Fig. O All the other ki ^ k2 blocks 
do not aff'ect the transfer efficiency rj (t) and average transfer 
time T at all. Especially, the initial component p^^,)t2 with ^ 
k2 will not influence 7/ (0 or r at any time t afterwards since 
it cannot evolve to the blocks with ki = fe. Therefore, only 
considering the dynamics of the (k, ^)-blocks are enough for 
the present purpose. 



V. TRANSFER EFFICIENCY AND AVERAGE TRANSFER 
TIME WITH CHANNEL DECOMPOSITION 

In this section we use the ^- space representation of master 
equation to calculate the average transfer time and transfer 
efficiency by the standard open quantum system method. As a 
highly organized array of chlorophyll molecules, the LH2 acts 
cooperatively to shuttle the energy of photons to elsewhere 
when sunlight shines on it. In this sense, we use the density 
matrix 



^^^^ " Z Z ^^^^'^'^2 ^^1 



(35) 



ki,k2 0,0'^A,B,C 



to describe the excitations in the initial state. From the discus- 
sions in the last section, only the ki = k2 = k blocks relevant 
to energy transfer. Therefore, there exists an equivalence class 
of initial states 

(0)] = {p I {O,k\p\O\k)=pok,O'k(0)} (36) 

that results in the same transfer efficiency and average transfer 
time as that forp (0). For further use, a special density matrix 
is chosen from the equivalence class ^(0) g (0)], 

^^^^ " Z Z ^^^'^'^ ^1 



k 0,0'^A,B,C 



(37) 



which satisfies {O,ki\g(0)\O\k2} = pok,0'k (0) for ki=k2 = 
k, and <0, ^i|^(0) fe) = for ki ^ ^2. ?(0) plays an 
equivalent role for determining the transfer efl&ciency and av- 
erage transfer time. Here, 



^^^(0)= poKO'kmo.k){()\k\ 



(38) 



0,0'^A,B,C 



is called as the ^-channel component of the density matrix. 
According to the above observation, we first choose every 
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(0) as the initial state to obtain the final state (0, which 
gives the /:-channel transfer efficiency at time t. 



ri''^ it) = r 



12 2 



(39) 



t O^B,C 

and the /^-channel average transfer time 



= ^ ^'Z Z ^o,,oAOdt\ (40) 



Then we prove a general proposition: 

For an arbitrary initial state p (0) (Eq. O LH2 com- 
plex, the transfer efficiency at time t and the average transfer 
time are the sum of ri^^^ (t) and r^^^ over all k-channels, re- 
spectively. 



;/(?) = J] ;/W (?) 

k 



T = 2_ ^ 

k 



[k] 



(41) 



In order to prove the above proposition we notice that the 
eff'ective initial state ^(0) evloves to 



(42) 



Since the corresponding transfer efficiency and average trans- 
fer time of ^(0) are 



k' 0=B,C 



Q0k',0t (f)df 



= I ^'Z Z ^ot^ot(^')dt\ (43) 

Eq. (HB is obtained from Eqs. gOl), gS]), and (l43l) . 

Namely, (0 and r are the sum of r]^^^ (t) and r^^^ for diff'erent 
momentum k channels. 

The present experimental observations^^ have provided 
some potential pathways for light-harvesting. One of them 
originates from the excitations on the B800 BChl ring. It 
shows that the excitations are transferred to the RC through 

B800 (LH2) ^ B850 (LH2) ^ B850 (another LH2) ^ > 

B875 (LHl) RC. As to our model, the initial state is spe- 
cialized as 

P (0) = PAkuAk, (0) |A, ki){A, k2\ . (44) 

kiM 

Accordingly, the ^-channel component of the eff'ective initial 
state ^(0) becomes 

(0) = PAKAk (0) |A, k) (A, k\ = pAkAk (0)?t^''^ (0) . (45) 

Taking 

^^'^^(0) = \A,k){A,k\ (46) 



as the initial state, we obtain the transfer efficiency and the 
average transfer time 

k' O^B,C 

^^'Xx,^t%(o^.^ (47) 

Hereafter, the superscript [A, k] denotes that the initial state is 
Eq. (l46l) . Similar to the above analysis about the proposition, 
we present a corollary: 

The transfer efficiency rj(t) and average transfer time r 
of the initial state in Eq. i44i are the weighted average of 
jj[A,.k] j-[A,k]^ respectively. 

ri(t) = YjPAk,Ak(0)r]^^^'^ (t) 

k 

r = YjPAk,Ak(0)T^^^'\ (48) 

k 

In the following, we will show the analytical and numerical 
results of t/"^^"^^ (0 and r^^'^^. 

First we consider the case without dephasing, i.e., F = 0. 
The time evolution from initial state g^^^^^ (0) only takes place 
in the (k, ^)-block. According to Eq. ([3Qb , the master equation 

of 'pk,k 



dpk,i 
dt 



= -i (HkPk,k - Pk,kHk) 



O^B,C 

gives the average transfer time 



- Yj {0l0kPk,k+Pk,k0l^k)^ (49) 



■^''' = = I^^'o^kiOdf- (50) 



When ^ = 0, Eq. (|49l ) about pok,0'k is rearranged as a sys- 
tem of diff'erential equations about vj (t) (j = 1 , • • • , 4) and 
V5 (0 = [V4(0]*: 



It is 



Vl = PAk,Ak^ 

V2 = PBk,Bk + PCk,Ck, V3 = PBk,Ck + PCk,Bk^ 
V4 = PAk,Bk + PAk,Ck, V5 = PBk,Ak + PCk,Ak- 



d 

— Vl (t) = ig+ [V4 (0 - V5 (01 , 

dt 
d 

-V2 (t) = -ig+ [V4 (0 - V5 (01 - (0 , 

dt 
d 

-V3 (0 = -ig+ [V4 (0 - V5 (01 - rv3 (0 , 

dt 
d 

— V4 (0 = (0 - ig+ [V2 (0 + V3 (01 

dt 



(51) 



2i (Ji - J2) + /AQ + 



V4 (0 , 



(52) 
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with initial conditions 

VI (0) = 1, V2 (0) = V3 (0) = V4 (0) = V5 (0) = 0. (53) 

Here g+ = (gi + 2g2). Solving the above differential equa- 
tions, we obtain 



n Jo 



g^ + (/i-/2+AQ/2)Hrg/4 



(54) 



with To = r/2, which is independent of the dimerization pa- 
rameter d. Similarly, when k = ±7i, the average transfer time 
of (^o) is 



^VA,k-.n^ ^ £ + (/i+/2^-AQ/2)^+rg/4 



(55) 



where g_ = {gi- 2g2). It is a quadratic function with respect 
to d. The optimal parameter d with the shortest transfer time 
satisfies 



^opt 



When ^1 = 2g2, Eq. (l55l) shows that r^"^'^"-^^ = oo, corre- 
sponds to 7/ = 0, the energy transfer is prevented at this time. 

If the dephasing is present, i.e., ^ 0, we can provide 
approximate solutions for r^"^'^"^^ and r^^'^"-^^, 

^[A.^=0L g+(4r.- r)/r+(27i- 272+ AQf+ 

_rA.^=..i_ g-(4r.- r)/r+(27i+ 272^- a^)^+ r^/4 

where = F -I- . They almost exactly agree with the numer- 
ical calculation below, and can also be confirmed by Eq. 
and (|55]) when F' = 0. The details are shown in Appendix C 



VI. ENERGY TRANSFER EFFICIENCY AND AVERAGE 
TRANSFER TIME IN NUMERICAL CALCULATION 

For a general k, the analytical solution of t/^^'^^ (0 and r^^'^^ 
is not easy to get. Nevertheless, the numerical results of t^^'^^ 
as a function of 6 are plotted as blue scatter lines in Fig. (H 
Here we have chosen 



h 
F 



0.3, 



1, 



^1 AQ F' 

^ =0.5, — =0.1,-- = 1, 

F F F 



(57) 



g2/F = for the upper panel, g^lY - 0.125 for the lower 
panel, and t is in the unit of (1 /F) and is long enough to ensure 
1] - 1 . It shows that when k ^ ^ and 6 varies from - 1 to 1 , 
there always exist optimum cases with and shorter 

average transfer time. This fact reflects the enhanced efl'ect of 
dimerization. 




FIG. 4. (color online) The average transfer time r^^'^^ of 5^^'^^ (0) 
(blue scatter lines) and Tmix of the initial mixed state p^.^(0) (red solid 
lines) with respect to the dimerization degree 6 of the B850 BChl 
ring. Here = 8, /i/F = 0.3, /2/r = 1, gi/F = 0.5, AQ/F = 0.1, 
FVF = 1, g2/F = (upper panel) and g2lT - 0.125 (lower panel). 
T is in the unit of (l/F) and rj - 1. It shows that each r^^'^^ {k ^ 0) 
curve has a minimum at (^jj^f ^ ^ 0. Tmix is the equal- weighted average 

of a complete set of jr^^'^^ }. 



We then take the mixed initial density matrix p(0) 
p^- (0) as an example. 



P:?>.x(0) = XipAM,(0)|A,;><A,y| 



= YjPAhAk2mA,h){AM- (58) 



The weight pAkM always satisfies 



1 ^ 1 



(59) 



From Eq. (l48t . the transfer efficiency and the average transfer 
time of p;*. is 



k 

T ■ - - V T^^'^^ 
Tmix - Zj ^ ' 



(60) 
(61) 



Tmix is also verified numerically and shown in Fig. |4]as the red 
solid lines. 

In order to see the dynamics of the transfer process clearly, 
we plot T/mix with respect to the dimerization degree d and time 
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150 



FIG. 5. (color online) (a) The contour map of the transfer efficiency 
of the initial mixed state //mixC^^, as a function of dimerization degree 
S and time t for the same setup as that in the lower panel of Fig. |4] 
(b) The profile of rj^^i^ along to = \2 (black dashed line in (a)), (c) 
The profile of rj^i^ along Sq = -0.5 (red dot line in (a)). It shows 
that r/j^ixi^, t) increases over time, and an optimum 6 can enhance the 
transfer efficiency. 




FIG. 6. (color online) Plots of r as a function of the dissipation ratio 
gi/r with J2IT = 0, 1,2,3, where /i/F = 0.3, 6 = (AQ/2 - /i)//2, 
g2lgi = 0.25, An/r = 0.1, rvr = 0.5, and ? is in the unit of (l/F). 
It shows that r decreases as gi/F increases, but increases with the 
increasing of /2/r. 



gives the contrary result. 

The minimal r^^'^^*^] 
(AD/2-/i)//2, 



is reachable at 



^opt 



t in Fig. Oa), i.e., //mix = ^mix(^, 0- At a certain instant = 12, 
^mixC^, ^0) as a function of 6 is plotted in Fig. [3b), while for a 
certain dimerization degree ^0 = -0.5, //mixC^o, as a function 
of t is plotted in Fig. Oc). Here the parameters are chosen as 
same as the ones in Fig. [4] except that g2lT = 0.125. The 
contour map Fig. [3a) and the profiles of y/mix(^, in Fig. [3b) 
and (c) show that (1) 7/mix(^, increases monotonously as time 
goes by. In the large t limit, y/mix(^, equals to 1. (2) At any 
certain short instant, an optimum S can enhance the transfer 
efl&ciency. 

Similar to //mix and Tmix, in general, there exists an optimal 
^opt ^ for an arbitrary initial p (0), which means that a suit- 
able distortion of the B850 ring is helpful for the excitation 
transfer. This result agrees with the x-ray observation that the 
Mg-Mg distance between neighboring B850 BChls is 9. 2 A 
within the ayS-heterodimer and 8. 9 A between the heterodimers 
reported in Ref.^^ The B 850 ring is indeed dimerized in na- 
ture. 

As shown in Fig. [H r^"^'^"^^ and j^^^^"-^^ are particu- 
lar since nearly all the other r^^'^^ are within the range of 
|^^[A,/:=o] .^[A,/:=±;r]j^ avcragc transfer time r of an 

arbitrary p (0). Besides, the absolute value of ^opf for 
the k = ±71 case is larger than the one of other p(0), i.e.. 



^opt| 



< s 



[A,k^±7T] I 

opt 



Hence, once we have known the proper- 



ties of T^^'^"^^ and r^^'^"-^^, the behavior of a general r can be 
conjectured to some extend. Compared the lower panel of Fig. 
[4] with the upper panel, a larger g2/T can increase r^"^'^"-^^ but 
decrease r^^^k"^\ in the g2/T = 0.125 case, the homogeneous 
pure state ^^'^"^^ (Q) is better than the mixed state p^.^ (0) for 
energy transport. However, the upper panel with g2/T = 



^[A,k^^n] ^ (gi-2g2f(4r + 3r) + r(r + r)^ /4 
2(gi-2g2)'r(r + p) 



(62) 



In the toy model illustrated in Fig. ^ J2 > Ji and g2 < gi. 
When 



we have r^^.'^^"^^ = r^^'^^^^ where 



(r + F) 



2 (2/2 - 2/1 - AQ) 



(63) 



(64) 



On the side of < g2/gi < 7g, r^^t'""^ < r^^'^^^^, while on 
the other side yg < g2/gi < 1, T^in^-""^ > r^^^^-^^ 

In general, the shortest average transfer time of an arbitrary 
initial p(0) is within the range of [T^'n""''^ '^^'^'^"^^]- The mean 
value of 7-^"^.'^"-^^ and r^"^'^"^^ can roughly reflect the influence 

rmn ^ 

of parameters on the transfer process. 



-hr' 



[A,k^O] 



)■ 



(65) 



In Fig. [6l we plot r with respect to gi/T for difl'erent J2/T = 
0,1,2,3. Here, Ji/T = 0.3, S = (AD/2 - /i) 7/2, gi/gi = 
0.25, AQ/r = 0.1, rvr = 0.5, and f is in the unit of (l/F). It 
shows that f decreases monotonously as /F increases. In the 
short gi/Y limit, r tends to infinity, which is reasonable since 
the two BChl rings are decoupled in this case. Moreover, r is 
larger when J2 /F is larger. 
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VII. CONCLUSION 



and its Hermitian conjugate 



In summary, we have studied the craggy transfer in Hght- 
harvesting complex with dimerization. We employed the open 
quantum system approach to show that the dimerization of 
the B850 BChl ring can enhance the transfer efficiency and 
shorten the average transfer time for different initial states 
with various quantum superposition properties. Actually our 
present investigation only focuses on a crucial stage in pho- 
tosynthesis - the energy transfer, which is carried by the co- 
herent excitations in the typical light-harvesting complex II 
(LH2). Here the LH2 is modeled as two coupled bacteri- 
ochlorophyll (BChl) rings. With this modeling, the ordinary 
photosynthesis is roughly described as three basic steps: 1) 
stimulate an excitation in LH2; 2) transfer it to another LH2 
or LHl; 3) the energy causes the chemical reaction that con- 
verts carbon dioxide into organic compounds. Namely, the 
excitations are transferred to the RC through B800 (LH2) 

B850 (LH2) ^ B850 (another LH2) ^ > B875 (LHl) ^ 

RC. Obviously, the first two are of physics, thus our present 
approach can be generalized to investigate these physical pro- 
cesses. Although photosynthesis happens in diff'erent fashions 
for diff'erent species, some features are always in common 
from the point of view of physics. For example, the photo- 
synthetic process always starts from the light absorbing and 
energy transfer. 

Another important issues of the photosynthesis physics 
concerns about the quantum natures of light^^'^^. Since the ex- 
periments have illustrated the role of the quantum coherence 
of collective excitations in LH complexes, it is quite natural 
to believe that the excitation coherence may be induced by the 
higher coherence of photon. Therefore, in a forthcoming pa- 
per we will report our systematical investigation on how the 
statistical properties of quantum light aff'ects the photosynthe- 
sis. 
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Appendix A: Equivalent non-Hermitian Hamiltonian 

In the case without dephasing, i.e., F^. = 0, an equivalent 
non-Hermitian Hamiltonian is introduced to study the dynam- 
ics of the open system, 

^ F- 

H = Hs-iY,i{B]Bj + C]Cj). (Al) 

M 

The equivalence between Eq. (lAll) and (fTTI) is shown as fol- 
lows. On the one hand, the Schrodinger equation 



Hs-iY,i{B]Bj^C]C,) 



7=1 



(A2) 



at 



Hs+i2i{B]Bj^C]C,) 



(A3) 



gives the evolution equation of the density matrix p = 



(A4) 



7=1 



On the other hand, when the dephasing terms are absent, 
the master equation Eq. (fTTI) becomes 



-£ = -i[Hs,p\ + Z r.lO^pOj - - {0]0j,p\]. (A5) 



7=1 O^B,C 



The above equation is written on the expanded Hilbert space 
with an additive vacuum basis |0). Compared with Eq. (IA4I) . 
the additive term in Eq. (IA5I) Zo=fi,c ^j^jP^j has only 
contribution to dpy^^/dt, which does not change the dynamics 
of the system. The Eqs. (IA4I) and (IA5b are equivalent for 
determining the time evolution of^oj,0'j'- 

For the non-Hermitian Hamiltonian, the corresponding 
transfer efficiency and the average transfer time are 

ri{t)= CfjTj |<0|O,|^(O)r^^' 

^0 j^l O^B,C 
^ poo 

T== ^^Tj^J Tj \<^\0M{n)\'df (A6) 

Due to the equivalence of the non-Hermitian Hamiltonian and 
the dissipative master equation, the results of Eq. (IA6I) are as 
same as the ones calculated by Eq. Eqs. ([T9l) and ([2Ql) . The 
non-Hermitian Hamiltonian method has an advantage over the 
master equation one for saving computer time. Instead of 
equations, only a system of N equations are needed to be 
solved in the non-Hermitian Hamiltonian case. 

However, when the dephasing terms are present, there 
is no equivalent non-Hermitian Hamiltonian. In this case, 
compared with Eq. (IA4I) , the additional term Yu0^b,c^) 
O^jOjpO^jOj cannot be omitted any more. It can also affect 
the evolution of the density matrix of the LH2 system. 



Appendix B: Transform the master equation to the ^-space 

In this section, we will transform the master equation from 
the real space (Eqs. O-Cl)) to the ^-space (Eqs. (I30l)-(l32l)). 
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Since Hg = Yjk and p is expressed as Eq. (1261) , we have 



- / [Hs,p] = -i Hk, ^ pk,,k2 + ^ (pv,k + Pk,v) 

k kiM k 

kiM 

-i ^ (i:^ytp/,,v - ^v,kfik) . (B 1 ) 



Here [Hs,pv,v\ = since they are in the different subspaces. 

According to the Fourier transformation Eq. ([22l) . the term 
2j OjpO\ becomes 



j kiMMM j 

= KMOk,Pkuk,Ol^ = Yj OkPKkOl (B2) 



kiM 



The term Y,j ^]^jP is transformed as 



kiMM j V ^2 



Pki,v 



- Y '^3 ^^1 + Pkl ,V 

kiM V k2 

= Y ^l^kiPkiM + Y ^l^kPk,v 



(B3) 



/^l,^2 

Similarly, 



kuk2 



Finally, the term Zj ^j^jP^j^j is written as 



•{k-i-k^+ks-ke)] q\ /o, , /ot 
i{h-ki+k2-k^)j r\\ /O, , /ot 



k2MM j 

= Y ^hM+k^-hOl^Ok.pk.MOlpk, 

kiMMM 

= — 2^ Oj^Ok,Pkuk20j^Pk2+k-ki' 



(B5) 



^1,^2,^ 



Therefore, Eqs. ([3Ql)-(l32l) are obtained by summarizing Eqs. 
(lBT1)-(lB5l). 



Appendix C: Approximative master equations for special 
cases 

In the cases of ^ = and k = ±n, we have the approximate 
master equation. 



dpk,i. 
dt 



= -i(HkPk,k-Pk,kHk)+ Y {^'OlOkPk,kOlOk 



O^B,C 



r + r 



{0l0kPk,k +Pk,kOl^k)\ ' 



(CI) 



It is verified numerically that the term 0\Okpk,kO\Ok in Eq. 
(ICll) plays the same role as (T' IN) Y,k' 0\,Okpk,kO\Ok' in Eq. 
^fork = 0, ±71. For the {k = Q,k = 0)-block, Eq. ([52) 
becomes 



Vl (0 = ig+ [V4 (0 - V5 (0] , 



V2 (0 



.[V4 (0-V5 (0]-rv2 (0, 



dt 
d 

dt 

^V3 (0 = -ig^ [V4 (0 - V5 (0] - (r + r) V3 (0 , 

d 

— V4 (0 = (0 - ig^ [V2 (0 + V3 (0] 



2/ (/i - /2) + + 



r + F 



V4 (0 , 



(C2) 



Solving the above diff'erential equation we have r^^'^"^^ shown 
in Eq. (1561) . The average transfer time r^^'^"-^^ for the k = ±7i 
channel is also obtained similarly. 
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